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Relaxation of a two-level system (TLS) into a resonant 
infinite-temperature reservoir with a Lorentzian spectrum is 
studied. The reservoir is described by a complex Gaussian- 
Markovian field coupled to the nondiagonal elements of the 
TLS Hamiltonian. The theory can be relevant for electromag- 
netic interactions in microwave high-Q cavities and muon spin 
depolarization. Analytical results are obtained for the strong- 
coupling regime, fio ^ i', where f2o is the rms coupling am- 
plitude (Rabi frequency) and v is the width of the reservoir 
spectrum. In this regime, the population difference and half of 
the initial coherence decay with two characteristic rates: the 
most part of the decay occurs at t ~ , the relaxation being 
reversible for t < (^o!/)"^^^, whereas for t > {Q.lv)''^^^ the 
relaxation becomes irreversible and is practically over. The 
other half of the coherence decays with the rate on the order 
of V, which may be slower by orders of magnitude than the 
time scale of the population relaxation. The above features 
are explained by the fact that sk, t <^ v^^ the reservoir tem- 
poral fiuctuations are effectively one-dimensional (adiabatic). 
Moreover, we identify the pointer basis, in which the reduc- 
tion of the state vector occurs. The pointer states are cor- 
related with the reservoir, being dependent on the reservoir 
phase. 

PACS numbers: 42.50.Ct, 05.40.-a, 76.75.-H, 03.65. Bz 



I. INTRODUCTION 

Decoherence of quantum systems is a subject of a sig- 
nificant current interest Q-J^, since this phenomenon is 
of a great importance both for the fundamentals of the 
quantum theory |^ and for the field of quantum compu- 
tation 1^ . In most treatments of the decoherence and dis- 
sipation the thermal reservoir have been assumed to have 
a sufficiently broad spectrum (a short correlation time), 
so that the system-reservoir coupling can be considered 
weak, which allows one to use standard approaches, such 
as master equations 

For the generic case of a two-level system (TLS) this 
approach yields the familiar relation between the popu- 
lation (or longitudinal) and coherence (or transverse) re- 
laxation times, Ti and T2 respectively: Ti > 12/2. The 
equality here is obtained in the case of a transverse reser- 
voir, i.e., when the reservoir variables enter only the non- 
diagonal elements of the TLS Hamiltonian. Otherwise, 
when the reservoir is coupled also to diagonal elements of 
the Hamiltonian, the decoherence may occur much faster 
than the population relaxation. 



Owing to the fast experimental progress in recent 
years, there appeared a number of novel materials and 
devices, such as photonic-band structures |^, high-Q cav- 
ities pO|-|l^, and semiconductor heterostructures, many 
of which are of potential relevance for quantum infor- 
mation processing. By creating structures in the elec- 
tromagnetic continuum, such devices often produce a 
strong system-reservoir coupling. In such cases the stan- 
dard master-equation technique is inapplicable and one 
should resort to other approaches. Fortunately, for zero- 
temperature reservoirs, the TLS relaxation can be cal- 
culated for the general case in quadratures jl^. How- 
ever, the finite-temperature case, which is important in 
the radio and microwave frequency ranges, is much less 
understood. Even the infinite-temperature limit, when 
the reservoir can be often modeled by a classical random 
field, is insufficiently studied. 

In this paper we report an analytical solution for the 
dynamics of a TLS which is strongly coupled to a trans- 
verse reservoir with a Lorentzian spectrum in the infinite- 
temperature limit. In this case the reservoir can be de- 
scribed by a classical field, which is a complex Gaussian- 
Markovian random process |l|,|l|. This problem can 
be of relevance for microwave high-Q cavities at 
temperatures of several degrees Kelvin and higher. It de- 
scribes also a relaxation of a two-level atom in a resonant 
chaotic field 

Moreover, the above problem attracted a significant 
interest recently |l^,|l^,^ in connection with NMR and 
muon spin depolarization experiments. In this case, nu- 
merical calculations showed that the squared coher- 
ence may decay slower than the population difference, 
violating thus the above inequality. However, the above 
numerical studies did not elucidate the dependence of 
the relaxation on the parameters of the problem, not to 
mention a physical explanation of the phenomenon. 

Below wc obtain a comprehensive picture of the TLS 
relaxation and provide a physical interpretation of it. In 
particular, we show that in the strong-coupling regime 
about half of the coherence still survives after the popu- 
lation relaxation is practically over! The results obtained 
here are applied to the discussion of the pointer states, a 
concept of a significant current interest in the theory of 
quantum measurements PP, ^ , ^^ . 

The paper is organized as follows. In Sec. || wc formu- 
late the problem and introduce an analogy between the 
Liouville equation for a TLS and the Schrodinger equa- 
tion for the spin 1. In Sec. |l| we discuss stochastic 



differential equations for partial averages. In Sec. IV 
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we present several limiting cases. The population and 
coherence relaxation is investigate d in Sees. [v| and Jvl , 



respectively, and discussed in Sec. VII. In Sec. |VII1 



obtain the pointer states and show that, in contrast to 
the previously reported results, they are correlated to the 
reservoir state. Section IX provides concluding remarks. 
The three Appendixes show details of the derivations. 



II. FORMULATION OF THE PROBLEM 

We consider a TLS with the states |1) and |2) and the 
resonant frequency Aq. The TLS Hamiltonian is 

H{t) = nAo\2){2\ ~ (h/2)[n,(t)\2){l\+R.c.]. (1) 

The interaction of the TLS with a high-temperature 
reservoir is described by a complex function Qdt) — 
u(t) + iv{t), which is supposed to be a complex Gaussian- 
Markovian random process with the correlation function 



(2) 



where fio is the rms coupling amplitude and is the 
correlation time of ^cit). The reservoir spectrum has a 
Lorentzian shape with the half width at half maximum 
equal to u. 

The Hamiltonian (|^) has been used to describe muon 
spin depolarization experiments |l9y2C| ]. This Hamilto- 
nian describes also [in the rotating wave approximation 
(RWA) and in the interaction representation] a TLS with 
the resonance frequency coupled to a Lorentzian reser- 
voir with the spectrum centered at uo^ (which is the case, 
e.g., for microwave high-Q cavities |l^Jll|]) or to a laser 
chaotic field with the nominal frequency lOc . The 

RWA validity conditions are f2oj l^ol ^ where now 
Ao = Wo — We- 
lt is convenient to write the TLS density matrix p{t) 
as the column vector, 



(3) 



where n = pn — p22 is the population difference and the 
superscript T denotes the transpose. Then the Liouville 
equation for p[t) can be written in the form 



where 



A{t) = i 



Ao -^l{t)/2 
-f7c(t) VLl{t) 
n,{t)/2 -Ao 



(4) 



(5) 



The solution of Eq. §) is r[t) = G[t)r{Q). Here G{t) is 
the Green function which obeys the equation 



G^A{t)G 



(6) 



with the initial condition G(0) — I, where / is the unity 
matrix. 



It is well known that the Liouville equation can be 
cast as the equation 123] 



s = s X B{t) 



(7) 



for the motion of a classical magnetic moment (pseu- 
dospin) 

s={s^,Sy,Sz) = (2Rep2i,2Imp2i,n) (8) 
in the effective magnetic field 



B{t) = {u{t),v{t), Ao). 



(9) 



This analogy helps to obtain an insight into the TLS 
behavior. 

Moreover, it is of interest to consider another anal- 
ogy, as follows. Under the linear transformation ip — Sr, 
where S" is a diagonal matrix with Su = —V2, Soo = 1, 
and S-i^-i = V2, Eq. (Q) becomes the Schrodinger equa- 
tion for a spin 1 in the magnetic field (|^), 



= iB{t) ■ Sip. 



(10) 



Here Si {i = x,y,z) is the operator of the ith compo- 
nent of spin 1 in the representation of the eigenfunctions 
of Sz (the cyclic-basis representation) |24|. [Equations 
(^) and (|o|) imply that the effective gyromagnetic ratio 
equals 1.] Note that the pseudospin s has the meaning 
of the polarization vector [Q for the above spin 1. The 
analogy ( |lO| ) can simplify calculations, by allowing one to 
use the standard textbook techniques developed for the 
Schrodinger equation. 



III. EQUATIONS FOR PARTIAL AVERAGES 

Consider the TLS density matrix averaged over such 
realizations of the random process ^dt) that assume the 
value 



ilc ^ u + iv = fie*'^ 



(11) 



at t. This partially averaged density matrix p{D,,t) writ- 
ten as r{(i,t) [see Eq. (||)], where (l — {u,v), is given 

by 



r{n,t) ^G{n,t)r{0). 



(12) 



Here G'(r2, t), the partial average of G{t), obeys the equa- 
tion lUl 



G ^ A{n)G + LG, 



(13) 



where A{fl) is given by Eq. (^ with a constant flc- The 
time dependence of the coupling is taken into account in 
Eq. (|lj) by the stochastic operator L = + Ly, where 
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The random vector Q, has a Gaussian stationary distri- 
bution dW{il.) = /(r2)c?ri, where = dudv and 



(15) 



Assuming that the TLS does not interact with the reser- 
voir at i < 0, the initial condition for Eq. ( |T3| ) is 
G(ri, 0) — If{^l), where / is the unity matrix. 
The fully averaged density matrix is given by 



where 



m = G(t)r(0), 



G{t) = / dnG{n,t). 



(16) 



(17) 



Alongside with the forward partial averages r(fl, t) and 
G(f2,i), it may be expedient to consider the backward 
partial averages f{Vl' ,t) and G{ft',t) related by 



fin',t)^Gin',t)r{o), 



(18) 



where the tilde denotes the average over (i{t) subject to 
the condition (l{0) — (l' . As shown in Appendix |^ the 
forward and backward partial averages are related by 

G{n,t) ^ s-^G^{n,-(f>,t)s^f{n). (i9) 

Henceforth we shall focus on the exact-resonance case, 
Ao = 0. (20) 
Consider a matrix Fourier series 



Here r„ is a diagonal matrix 
r„ = diag(e 



i{n+l)4> ^in4> gi(n-l)0-j 



(21) 



(22) 



and G„(r2, t) — (T„G(f2, i))^, where the average over (j) is 
denoted by (. . .)^ = (Stt)"! dep... . Muhiplying Eq. 
( p^ ) by T„ from the left and integrating the both sides 
of the resulting equation with respect to from to 27r, 
one obtains the equations 



G„ = A{n)Gn + M„G„ 



(23) 



where A{n) = ^(17,0 = 0) and M„ = 
diag(L„+i, L„, i„_i). Here Lk is defined so that for an 
arbitrary function F{n) 



(24) 



On writing L in the polar coordinates p9|] , one obtains 
that 



in - 



2 an^ 



'2n 



(25a) 
(25b) 



Taking into account that the only nonzero initial con- 
ditions for Eqs. (H) are G„;^_„(f], 0) = f{n) {n = 
0, ±1), one obtains from Eqs. (P3| ) and (|2l| ) that 



G{n,t) 

R{fl,t) 



e-">'E{n,t)/2 e-^^'l'P{n,t) 



e^^'t'P{n,t) -e"l'E{n,t)/2 R{n,t) 

(26) 

The functions entering Eq. (|2^) satisfy the following sets 
of equations, 



N^-inE + LoN, E^-iVLN + LiE 



(27) 



and 



P = -{i/2)nQ + L2P, Q = in{R-P)+LiQ, (28a) 
{i/2)nQ + LqR. (28b) 

The only nonvanishing initial conditions for Eqs. (|2 
and (|2l) are 



(a) N{n,o) = f{n), (h) R{n,o) ^ f{n), 



(29) 



respectively. The functions N{n,t), R{n,t) and P{fl,t) 
are real, whereas E{n,t) and Q{D,,t) are purely imagi- 
nary. This follows from the fact that, on changing the 
variables E'{n,t)_=^ iE{Vl,t) and Q'{^,t) = iQ{n,t), 
Eqs. ( p7|) and (|2^) become sets of equations with real 
coefficients and real initial conditions. 



(30) 



Averaging Eq. (12) over yields 
r{n,t)^G{rt,t)r{0), 



where r{n,t) = {r{Q,,t))cj,. As follows from Eq. (12£ 

G{n,t) = {G{n,t))^ = dmg{R{n,t),N{n,t),R{n,t)). 

This means that the population relaxation and the co- 
herence relaxation proceed independently of each other. 



n{n,t) = N{n,t)n{0), 



(31) 
(32) 



Hence, finally, 

nit) = N{t)n{0), pi2(2i)(t) = i?(0Pi2(2i)(0), (33) 

where the population and coherence relaxation functions 
are obtained by 



N{t) = 2tt / N{n,t)ndn, 

Jo 

/>oo 

R{t) ^2tt R{n,t)ndn. 

Jo 



(34a) 
(34b) 



3 



IV. LIMITING CASES 



(39) 



a. Zero temperature. For the sake of comparison, 
we niention here the results obtained for the zero- 
temperature reservoir with a Lorentzian spectrum [^ . 
In this case, for an arbitrary initial state the TLS relax- 
ation is determined by the equalities P22{i) = No{t)p22{Q) 
and pi2{t) = Ro{t)pi2{Q). The functions N^{t) and Ra{t) 
can decay monotonously or in an oscillatory fashion, de- 

f ding on the coupling strength, as described in Ref. 
. However, for all values of the coupling strength and 
detuning one can show that A''o(t) = |-Ro(i)P, which is the 
extension of the relation Ti = T2/2 for a nonexponential 
relaxation. The above equality is in a sharp contrast with 
the present infinite-temperature shown below. 

6. Weak coupling. Returning to the infinite- 
temperature case, two limiting cases can be readily con- 
sidered. For a weak coupling, f^o ^ the average popu- 
lation difference and coherence decay exponentially with 
the decay times Ti = v/ilf^ and T2 = 2Ti respectively 



c. The static limit. In the opposite case v ^ (th e 
static lirnit) one can set L„ « in Eqs. ( p7| ) and (28), 
yielding ||,0,||,|6| 



N,,{t) ^i-n„tF{not/2) 

Rstit) = [iVst(i) + l]/2, 



where F{z) is Dawson's integral 



F{z) 



dye^ 



(35a) 
(35b) 



(36) 



The function (|35a| ) (see Fig. |l|, the dot-dashed line) sim- 
plifies in two limits. 



N,t{t)^i-nlty2 {t^n^'), 
N,t{t)^~2/nie {t:$>n^'). 



(37a) 
(37b) 



It vanishes a.t t = 1.85f2Q and has one minimum equal 
to -0.285 at t = 3.001^0 . 

As follows from Eq. ( |35bD , limi^oo Rstit) = 0.5. This 
nonzero limit results from the fact that in the static limit 
the component of the TLS pseudospin parallel to the ef- 
fective field B = {u,v,0) is conserved. Equation ( |35| ) 
implies that for t <C one has N(t) ^ R'^it), as in the 
weak-coupling case. In contrast, for t 3> ^Iq^ the ratio 
\N{t)\/R'^{t) tends to zero. 

Note that in the static limit, the relaxation occurs due 
to the statistical spread of the coupling amplitude and, 
correspondingly, in principle it can be reversed by an echo 
technique. Irreversible relaxation is obtained only in the 
presence of temporal fluctuations ^ 0). 

d. Short times. Consider effects of the temporal fluc- 
tuations in the strong-coupling regime, 

rJo > 1^, (38) 
for sufficiently short times, i.e., as shown below, for 



Here D = is the diffusion coefficient for Q,{t) 

[cf. Eq. (0)]. For short times, one can use the time- 
dependent perturbation theory, as described in Appendix 
to calculate the backward partial average G(5l,t), 
which, in turn, yields G{Q.^t) by Eq. (p^. Comparing 
the resulting expression with Eq. (E^) yields 



N{Vl,i) w f{Vl){a{n,t) coaftt 

+ [D{1 - V?t^)l{2V?)\s\i^m}, 



(40) 



R{n,t) 



N{n,t)/2 

+ cos ilt) - 



^ f{n)[i/2- {Dt/2n'^){2 

{3D/2n^)smnt], (41) 



where a{n,t) = 1 - Dt^/S - Dt/2n'^. Expressions for 
E{il,t), P{n,t), and Q{il,t) are shown in Appendix ^ 



V. POPULATION RELAXATION 

First, we discuss the population relaxation. Equations 
( p7| ) are related to the Schrddinger equation for some 
TLS {\a),\b)} coupled to a resonant chaotic field, as fol- 
lows. Writing the Hamiltonian of the latter TLS in the 
form 



Ho = hVc(t)\b)(a\+R.c. 



(42) 



where Vc(t) is a complex Gaussian-Markovian process, 
one obtains the Schrodinger equation 



Ua 



-iV*{t)uba, Uba = ~iVc{t)Ua 



(43) 



where Uaa{t) and Ubait) are matrix elements of the evo- 
lution operator U(t). If the Hamiltonians (|^) and ( |4^ ) 
describe the same interaction, then Vc{t) = —flc{t) /2 and 
hence 



v = n/2, Vo = no/2, 



(44) 



where V = \ Vc\ and = (V^). 

The averaging of the above quantity U{t) was discussed 
in Ref. [^. Writing equations for the partial average 
quantities UaaiV,t) and Uba{V,t) ^ Eqs. (3.3), (3.7)] 
and eliminating the field phase, one obtains the equations 
pO| which can be identified with Eqs. (|2^) under the 
substitutions 



Uaa{V,t) = Uaa{V,t) ^ N{n,t), 

Uba{V,t)e''^ E{n,t), 

(a) v^n, (b) Vo ^ no. 



(45) 
(46) 
(47) 



For instance, for <C iloi Eq. ( ^0| ) can be obtained from 
Ref. H, Eq. (5.25)], in view of Eqs. (H), (||), and dH). 

As follows from Eq. (45), the quantity Uaa{t) ~ 
JdVUaa{V,t) equals N{t) under the substitution (|^). 
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Thus, one can apply directly the results of the compre- 
hensive study of Uaa{t), performed in Ref. to the 
population relaxation function N{t). Plots of N{t) for 
different values of i^/Qq are shown in Fig. |l]. 
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FIG. 1. The population relaxation function N{t) versus the 
dimensionless time Qot for the values of v/^lo shown in the 
plot. 

Henceforth, we focus on the strong-coupling regime, 
Eq. (^8|). In this case, as follows from , the population 
relaxation is described by 



N{t) ^ N,t{t)J{at), 



(48) 



where a = (2£>)^/^ — (iloJ^)^/-^ and J{at) is a dimen- 
sionless function of at, which describes the irreversible 
relaxation. For short times, the function N{t) is pit] 



J(at) w 1 - Dt^/3. 



(49) 



One can show that the second-order correction to Eq. 
( p9| ) is on the order of D^t^, which implies the validity 
condition (^8|) for the above short-time results. 

For 3> D^^ J(at_) tends to zero, performing damped 
oscillations (see Fig. g). Correspondingly, the character- 
istic rate of the irreversible relaxation is on the order of 



VI. COHERENCE RELAXATION 

Consider the coherence relaxation. We shall focus on 
the strong- interaction regime, Eq. (p8|). At short times 
the coherence relaxation function R{t) is obtained by in- 
serting Eq. (41) into ( t34b| ) and performing the integra- 
tion. For t ^f^o ^ the R{t) is very close to the static 
result, the discrepancy increasing with t. As discussed in 



Appendix ||, for Q^'^ < < D^^/s 

R{t) « 1/2 - l/{nlt^) + vt{Co - Innot). (50) 
Here Cq — 5/3—7/2 = 1.38, where 7 is the Euler constant 




at 

FIG. 2. The function J(at) versus the scaled time at. 
The coherence relaxation for t^ ^ D^^ can be found. 



as follows. Casting Eqs. (|28|) as an equation for the 
column vector 



(51) 



and performing the linear transformation x(rj,t) = 
Sq{^l,t), Eqs. (ESh become 



X = [ins^ + M)x. 



(52) 



Here M is a diagonal matrix with Mmm = ^i+m- Note 
that Eq. (^2|) can be obtained also from partially aver- 
aged Eq. (^, ^{n,t) = iB-Si} + L^I;,withB = (m,w,0), 
on defining 



Xm{n,t) = (e^(i+")'^^„(j1,i))^/V^_i(0) 



(53) 



(m = 1,0,-1). Thus, x(^j^) has the meaning of a 
weighted partial average of the wave function for spin 
1 in the stochastic magnetic field B{t) = {u{t),v{t),0). 

Next, we invoke the semiclassical dressed-state repre- 
sentation, by diagonalising the first term on the r.h.s. 
of Eq. (p^). Since Eq. ( ^ has the form of a (par- 
tially averaged) Schrodingcr equation, this can be done 
with the help of the standard finite-rotation operator. 
Namely, we rotate the z axis around the y axis by 7r/2, 
Xlf^,t') = Sox{^,t), where So = is defined in 

[02i. As a result, one obtains 



X' = [ins, - {D/2n^)C + Lo]x'. (54) 



Here 



3 -\/2 1 \ 
C = I -V2 4 -x/2 . (55) 
1 -%/2 3 / 

The vector x' is related to qhy x' — Siq, where 

-l/\/2 I/V2 I/V2' 
Si^SoS= \ 1 1 I . (56) 

-I/V2 I/V2, 
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Correspondingly, the initial condition for Eq. (jsl) is 
[of. Eq. (lib)] x'i^,0) = fmx^°\ where = 
(l/V2,l,l/\/2)^. 

The inverse linear transform q = Si ^x' yields 



R{n,t)^Rex'i{n,t)/V2 + x'o{^,t)/2. 



(57) 



Here we took into account that y\ (il,t) — x'*-i{^jt)^ 
which follows from the form of Eq. (jsj) and the initial 
condition for it. On introducing new variables by 



xLi^,t) ^ K^in,t)x,'^\ 

Eq. ( |57| ) becomes 

R{n,t) = [ReKi{n,t) + Ko{n,t)]/2. 



(58) 



Note that K_i{n,t) = K^{n,t). As follows from Eq. 



(60) 



For 3> g, the nondiagonal terms of the matrix coeffi- 
cient in Eq. (|5J) are much less than the differences of the 
diagonal terms (which are on the order of 2fl) and can 
be neglected in the first (secular) approximation. This 
results, in view of Eqs. (pq), in the following equations. 



ki = {in - 3D/2n^ + La)Ki, 
ko^i-2D/n^ + Lo)Ko, 



(61a) 
(61b) 



with the initial condi tions (|60|). 

Co nsid er now Eq. (61a). Note first that in Eqs. (61a) 
and ( |25b| ) the second terms in the parentheses are small 
as compared to the first terms, respectively, and hence 
can be neglected. Furthermore, performi ng t he change 
of variables Ki{Q,t) = e''^K{Q,t) Eq. (61a) becomes 
approximately 



k^{in + Ln)K, 



(62) 



where Lq is given by Eq. (|lj) with u ^ 17. Equation 
( |6^ ) describes the dephasing of a two-level system due 
to Gaussian-Mar kovian frequency fluctuations. The solu- 
tion of Eq. ( |6^ p3| implies that in the strong-interaction 
regime ( |38| ) considered 

Ki{n,t) « K{n,t) « /(rj)e'"*-^*'/3, (63) 

i.e., Ki{n,t) decays with the rate on the order of a. 

On the other hand, for < a, the equations (^4|) for 
Xm(ll,t) are strongly coupled. An analysis of Eq. ( |5^ ) 
shows that for SI < a, x'(^i^) ^^nd hence K„i{fl,t) [see 
Eq. dSS])] disappear on the time scale a~^. As a result, 
in view of Eqs. (||) and (|6|), for t > 



0, n<a 

Ko{n,t)/2 17 > a. 



(64) 



Since a <C r^o, the fact that KQ{n,t) « for 17 ^ a can 
be taken into account approximately by the boundary 
condition to Eq. (61b), 



(65) 



Equation (61t) can be solved by the conjecture 

Ko{n,t) ^ f{n)g{h{t)n), (66) 

where the functions g{X) and h{t) are to be found, using 
the initial and boundary conditions ( |60| ) and (65). As 
shown in Appendix ^ the solution yields 



Ko{n, t) - C2/(17)C'A/(fc, 2k + 1, -C), 



(67) 



(59) where k ^ 1/V2, C = n? inUe"^"^ - 1), M() is the degen- 



erate hypergeometric function [p7| , and 

C2 = r(i + fc)/r(i + 2fc), 



where r() is the F-function 
As follows from Eqs 




where 



4H) and (B4), 



R{t)^Ko{t)l2 (t»a-i). 



POO 

Km{t) ^2tt K^{n,t)ndn. 

Jo 



(68) 



(69) 



(70) 



Inserting here Eq. ( pTI) and performing the integration 
with the help of Ref. ||3J, Eq. 7.621.4], one obtains 



Ka{t) = Cie-'""''F{k, k; 1 + 2k; e"^'"*) 



(71) 



Here Ci = r2(l -|- fc)/r(l + 2k) « 0.66 and F()is 
the hypergeometric function The function ( [7l| ) 

monotonously decreases from 1 to with the average 
rate on the order of v. Equations ( |69| ) and ( [n] ) yield the 
following limits 



R{t)K 1/2 +{iyt/2){\ni^t-C3) (c 



R{t) « Cie 



-2kut 



(e^-^* » 1). 



<i<zy"^), (72a) 
(72b) 



Here C3 = 1 + 2k - 2^p{l + k) - 2j - ln2 « 0.14, where 
■00 is defined in Equation (|72^)was obtained with 
the help of the expansion 15.3.11 in p7| for the hyperge- 
ometric function. 

Equations ( |50| ) and (72a) describe the behavior of R(t) 
at t ^ The last term on the right-hand side (rhs) of 
Eq. ( |50| ) (which arises due to the temporal fluctuations) 
decreases with the time, increasing by the magnitude, 
in contrast to the second term (pertaining to the static 
limit). Note, however, that the last term on the rhs of Eq. 
( pO| ) is much less than the second term, i.e., R(t) ~ Rstit) 
for t <C a~^. The magnitudes of the above terms become 
of the same order (i//r2o)^^^ for t ~ , where R{t) has 
a maximum. Finally, for t ^ decreases, the 

coherence dynamics being determined by the temporal 
fluctuations [see Eqs. (||), (0), (^]. 
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In the first approximation, R{t) can be described for 
all times by the formula 



R{t) ^ R,t{t)Ko{t). 



(73) 



Indeed, as follows from Eq. (72a) and the above discus- 
sion, for t ^ the error in Eq. ( [73| ) increases with 
t from to the value of the order of (v/flo)'^^'^ ^ 1 at 
t ~ a~^, whereas for t Eq. ( fz^ ) is close to ( |69| ) 

with the relative error not exceeding by the order of mag- 
nitude {i^/nof^^ < 1- 

Note that the interpolation formula ( |73| ) is not unique. 
Instead of it, one can use with the same accuracy, e.g., 
the formula 



Kit) = [Nst{t) + Koit)]/2. 
The both formulas provide similar results. 



(74) 



DC 




1 2 

FIG. 3. The coherence relaxation function R{t) ver- 
sus logjQ f2o^ for the values of u/Qo shown in the fig ure. 
Solid lines, numerical solution; dashed lines, Eq. ([73[) for 
u/Qo = 0.01 and 0.0001. 

Figure H demonstrates R(t) for different values of ^o/^- 
The solid lines are calculated numerically by inverting the 
continued fraction describing a Fourier transform of R{t) 
pof , whereas the dashed lines are the plots of Eq. ( |7^ ) 
for v/flo = 0.01 and 0.0001. As follows from Fig. |, Eq. 
( [73| ) provides a good approximation to the exact result in 
the strong-coupling regime, the accuracy increasing with 
the decrease of ly/ilo (in particular, at v/ilo — 0.0001 
the solid and dashed curves coincide almost completely 
in Fig. I). 



VII. DISCUSSION 

As shown above, in the regime of a strong coupling to 
the reservoir, fio ^ J^, the TLS behavior is rather compli- 
cated. The population relaxation is characterized by two 
time scales. On the time scale ^q^, where most of the 
population relaxation occurs, the latter is of a reversible 
character, whereas on the time scale a^^ the population 



relaxation becomes irreversible and proceeds to comple- 
tion. 

The coherence relaxation is even more complicated. 
Roughly speaking, half of the coherence decays similarly 
to the population relaxation. However, for t ^ a~^, 
when the population relaxation is already completed, al- 
most half of the initial coherence still survives. The de- 
cay time of the latter is on the order of h'~^. This time is 
greater than the largest population relaxation time 
by the large factor {flo/i^)'^^^. 

To explain the above behavior, we write the pseudospin 

as s = sj_+S||, where sj_ = {Sx, Sy, n) and sy = (si, s|, 0) 
are respectively perpendicular and parallel to the initial 
effective field B{0) = («', v' , 0). Here we denoted ^^0) = 



= n' 



J4>' 



The pseudospin component s± 



describes the population difference n and the "out-of- 
phase' coherence = {s^ +iSy)/2 = ie^' 



'Im(p2ie-^'^'), 



whereas S|| (t) describes the "in-phase" coherence P21 — 

{si + isl)/2 = e"f''Re{p2ie-'^'). 

Except for the case il' ^ a, which has a negligibly 
small probability, during the time ~ the direction 
of ri(t) almost does not change. As a result, the back- 
ward average of s±{t) (i.e., the average with a given fi^) 
and hence the backward-averaged population and out- 
of-phase coherence decay at the rate a, as in the case of 
coUinear (adiabatic) field fluctuations l7|, |29| , |33|| . By the 
same reason, the backward-averaged in-phase coherence 
almost does not decay at t < a~^. Actually, it can decay 
significantly only after the field rotates by an angle of the 
order of 7r/2, which, on the average, requires the time of 
the order of the correlation time h'~^. The averaging of 
the above quantities over fi^ (i.e., the full averaging) does 
not change the time scales, at which they vanish, whereas 
it provides equal contributions from the in- and out-of- 
phase coherences. This explains the above fact that the 
populations and half of the coherence relax faster than 
the other half of the coherence. 

The above argument can be summarized as follows. 
The reservoir considered here is two-dimensional and 
symmetric with respect to the rotations in the xy plane 
(of the pseudospin space). However, for t ^ v'^ it be- 
haves effectively as a one-dimensional one, since the di- 
rection of the effective field B{t) does not change signif- 
icantly over short time intervals. In the strong-coupling 
regime, this fact results in a substantially asymmetric 
behavior of the pseudospin: the component s±{t) decays 
much faster than s\\{t), as discussed above. In contrast, 
in the weak-coupling regime, where the decay rates Ti 
and T2 are much longer than the correlation time 
(see Sec. IV), the short-time behavior of the reservoir 
is not important, and the density matrix elements relax 
with similar rates. 
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VIII. POINTER STATES 

The above results have a bearing on the the pointer 
states and related concepts [^P pl| , p^ . It is well known 
(Ref. and references therein) that an important part 
of any quantum measurement is the reduction of the state 
vector, i.e., the diagonalization of the density matrix of 
the measuring apparatus in a special (pointer) basis, the 
diagonal elements remaining intact. The apparent con- 
tradiction of this process to quantum mechanics, where 
any evolution is unitary, has not obtained yet a com- 
pletely satisfactory explanation. The resolution of the 
above contradiction suggested by Zurek |^ is that the 
reduction of the state vector occurs due to the coupling 
of the apparatus with the environment. 

One can ask which conditions the coupling of the ap- 
paratus (called below the system) and the environment 
(reservoir) should satisfy to be capable to produce the 
state reduction. There are, at least, two such conditions, 
as follows. The first condition requires that the pointer 
basis be independent of the state of the reservoir pi] . 
Another condition stems from the apparent contradic- 
tion between the role, which the reservoir, according to 
Zurek, plays in obtaining the state reduction, and the sec- 
ond law of thermodynamics, which postulates that any 
interaction with the environment ultimately results in the 
thermal equilibrium of the system, where, of course, no 
information of its initial state is preserved. The above 
contradiction can be removed by requiring that in the 
pointer basis the relaxation of the off-diagonal elements 
of the density matrix proceed much faster than the re- 
laxation of the diagonal elements. Then there exists a 
nonvanishing time interval, during which a measurement 
can be performed. Hence the second condition states 
that the system-reservoir coupling should be such that 
the relaxation of the system to the equilibrium proceed 
with significantly different time scales. 

Until now, the existence of the pointer states has been 
verified for the harmonic oscillator and a collection 
of many two-level systems Q weakly interacting with the 
environment. As concerns a single TLS, a fast phase re- 
laxation (T2 <C Ti) is known to produce the state reduc- 
tion in the basis of the energy eigenstates, as was noted 
l^sl in connection with the quantum Zeno effect 

As shown above, in the strong-coupling regime ( p8^ the 
relaxation of a TLS proceeds with significantly different 
time scales, which suggests a possibility of the e xist ence 
of the pointer states. 



VII im 



The discussion in Sec. 
plies that likely candidates for the pointer states are the 
eigenstates \ip±) of the initial value of the Hamiltonian, 
-{h/2)[fl'^\2)(l\ +H.C.]. One can easily check that 



|V;±) = (|l)±e^^'|2))/V2. 



(75) 



To prove that the states (75) are indeed the pointer 
states, one should consider the average density matrix 
conditioned to a fixed value of the initial field phase (j)' , 



p{4>\t)^2TT p{n',t)f{n')n'dn'. (76) 

Jo 

Consider the vector 

f(^',t)=G(0',i)r(O). (77) 

The function G(</>',<) = 2Tr G{n' ,t)f{n')n'dn' is re- 
lated to 



G{<j),t) = 2n / G{n,t)ndn 

Jo 

by [cf. Eq. (|l|)] 

Combining Eqs. (|6|), ^ and (^) yields 



(78) 



(79) 



G(0',i) 



R{t) -e-"f''Q{t)/2 e- 



N{t) 



e2"^P(t) e"^Q(t)/2 



'm 

R{t) J 



-i<p' 



(80) 



The functions entering Eq. (30) are obtained from the 
functions appearing in Eq. (26) by the integration, as 
in Eqs. (|3^). One can express q{t) = (Plt),Q{t), R{t)Y 
[cf. Eq. (|l|)] through K„,{t), Eq. (|7^), using q{t) = 

S^\'{t), where x'„M ^ K,n{t)x'm^ [cf. Eq. ^)]. This 
yields 

P{t) = [2Ko{t) - Ki{t) - if-i(i)]/4, (81a) 

Q{t) = [i^i(i) - K_i{t)]/V2, (81b) 

Kit) = [2Ko{t) + Ki{t) + i^-i(t)]/4. (81c) 

Inserting Eq. ( ^0|) into ( [77[ ) yields f{4>',t) and hence 
p{<j)',t) [cf. Eq. ^Spj. The transformation of the latter 
to the basis (^ is performed by p'{<j)',t) — p{<j)' ,t)S , 
where 



V2\e'* 



1 



(82) 



Finally, one obtains the following expressions for the com- 
ponents py((/)',i) +,-) of p' {(/)', t), 

p++(__)(0', t) = 1/2 ± Ko{t)Re[p2i{0)e~"^'], (83a) 
p+_(0', t) = p%{4>', t) - [N{t) + Q{t)]n{Q)l2 
+ilm{[Ki{t) + K^i{t) - 2£;(t)]p2i(0)e~''^'}/2. (83b) 

Equations ( ^3|) imply that p' {((>', t) becomes diagonal 
at t 3> since the functions of time appearing in rhs 
of Eq. ( 83b ) vanish at such times [see Sees. 0and0 
and Eq. ( |81b| )]. Taking into account that Ko{t) ~ 1 for 
t , in the time interval 



(84) 
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one obtains that 



P++(0) 
p__(0) 



(85) 



where )(0) = 1/2 ± Re[p2i(0)e"*'^'] are the diag- 

onal elements of /5' (</>', 0). 

Thus we have proved that in the interval ( p4[ ) the reser- 
voir produces the reduction of the state vector in the basis 
([75|). In contrast to the previous cases described in the 
literature, here the pointer states ( |75| ) depend on the state 
of the reservoir, due to their dependence on the reservoir 
phase (j)' . Since this violates the above first condition, the 
strong coupling considered here is not of a type allowed 
for the apparatus-environment interaction. 



IX. CONCLUSION 

Above we presented a comprehensive analysis of the 
relaxation of a TLS coupled to an infinite-temperature 
reservoir. We showed that in the strong-coupling regime 
the decoherence can proceed much slower than the pop- 
ulation relaxation and provided a physical interpretation 
of this effect. We identified the pointer states, which, in 
contrast to the previous findings, appeared to be corre- 
lated with the reservoir. 

The present results can be checked, e.g., in experiments 
with high-Q microwave cavities [^,^. Since the theory 
holds in the infinite-temperature limit, one should require 
that, at least, the average number of the photons Tiph in 
the resonance mode be large. Note that for a cavity mode 
with the frequency 21.5 GHz, used in experiments in Ref. 
I p^ , at the temperature 5 K riph — 4,A. An increase in 
the cavity temperature T and/or mode wavelength Ac can 
significantly increase the above number, since for high 
temperatures (nph S> 1) nph TAc. 

A more detailed estimation of the experimental condi- 
tions would require a consideration of corrections to the 
above results due to finite, though large, values of nph, 
which is out of the scope of the present paper . We 
believe, however, that the main results of this paper, in 
particular, those concerning the anomalously slow deco- 
herence and the pointer states, will remain valid, at least, 
qualitatively, also for moderately large Uph- 
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APPENDIX A: RELATION BETWEEN 
FORWARD AND BACKWARD PARTIAL 
AVERAGES 

Let us derive the relation between G{Cl', t) and G(J7, t). 
It is convenient to consider first the Green function U{t) 



of Eq. (0) defined by ^l){t) = U{t)^{0). The back- 
ward and forward partial averages of U (t) are denoted 
by U{(l',t) and U{(l,t). The latter function satisfies the 
equation 



il = [Ai{n)+L]U 



(Al) 



with the initial condition U{ft,0) = If{fl), where 
Ai{D,) — iB ■ S. Equations for backward averages were 
derived in Appendix A of Ref. p^ . In particular, Eq. 
(All) in Ref. ||2^ implies that C/(j1,t) obeys the equa- 
tion 



dt 



[Al{n) + L^]U^ 



(A2) 



with C/(f7,0) = /. 

The Gaussian-Markovian random process satisfies the 
detailed-balance condition, which implies that Eq. (A2) 
can be recast as 



dt 



(A3) 



where f/i(f7,t) = lf{n,t)f{n). Note also that Ai{n) = 
Ai{n, (j)), as defined above, obeys A'[{U, (f) — Ai{Q., — 



With the account of this, the comparison of Eqs. (Al) 
and (A3) shows that 



u{Q.,t) = u^{n,-4),t)f{n). 



(A4) 



Taking into account that U{n,t) = SG{n,t)S^'^ and 
U{n,t) = SGin,t)S-^, Eq. (^) yields finally Eq. (|l|). 



APPENDIX B: PERTURBATION THEORY 

Consider the perturbation theory in the strong- 
coupling regime, Eq. (|38|). One can write fl(t) — 
n' + W{t), where (l' = {u',v') = n{0) and W{t) ^ 
[W^i(0: W^2(0] obeys H^(0) = 0. Under the unitary trans- 
formation V'' = S'2'0, where Si^ram' = D^^,(0',7r/2,O) = 
d^,(7r/2)e™''^' is the finite-rotation matrix |3|], Eq. 
( p^ ) becomes 



^' = i[n's, + A2{t)]4>'. 



(Bl) 



Here A2{t) = -~2[W[{t)S^ + VF^(t)5j,], where W[{t) 
and Wl^it) are defined by W[{t) -f iWl^it) = [Wi{t) 

lW2{t)]t 



-i<f)' 



Solving Eq. (Bl) with the help of the time-dependent 
perturbation theory of the second-order in A2 (t) and av- 
eraging the result over W'{t) yields 



ip'in^t) = u'{n',t)ii;'{o), 



(B2) 



where 
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Jo 

ft l-tl 
/o Jo 

xU^{t2). (B3) 

Here the tilde and the angular brackets denote the aver- 
age over W'{t), whereas Uq^I) — e*^'*'^^. 

The functions 17' + Wl{t) and W2{t) are independent 
identical Gaussian-Markovian processes, which equal re- 
spectively fl' and at < = 0. For t <C W^2(0 can 
be considered as a diffusion (Wiener) random process. 
The same holds for W + Wl{t), unless fl' is too large 
in comparison with ^Iq. As follows from Eq. (py), the 
approximation (B3) holds for t <C a^^ ^ which 
means that for fl' < fig W{{t) and (i) in Eq. (^) 
can be considered as independent identical diffusion pro- 
cesses with the initial conditions ]¥[ 2(0) = 0. Perform- 
ing the averaging in Eq. (B3) with the help of the rela- 
tions {Wlih)) = and {Wl{ti)W^{t2)) = 2Dt26^j, where 
6ij is the Kronecker symbol, yields 

u'{n',t) = u^{t)[i - {Dt^/s)sl] -2D f dh r dt2t2 

Jo Jo 



xu^{t-h)s.Mh-t2)SyU'^{t2). 



(B4) 



Under the inverse linear transformation, f(J7',t) 



S^^^'{n',t), where 5*3 = S2S, Eq. (|2|) becomes (|l 



with 



G{Q.,t) = S^^U'{Q.,t)S^. 



(B5) 



From Eqs. (B4), (B5), and (p^ ) one obtains finally 
G{n,t) and hence Eqs. ®, QT^and 



E{n,t) w -2if{n) 



a(n,t) . ^ Dt^ 
sm S It H cos i It 

2 4n 



, (B6) 



+ 



1 1 
i:>(i - 



1 - 



Dt' 



3Dt\ 



■ sin fit 



2172 y 

Dt 
172 



cos 17t 



(B7) 



Q(17,i)«z/(17) 



L>i3 SUA 
" ~ " 2172 j 
i:»(172t2 _ 4) 



iflt 



2173 



2L> 

cos S 2i -I- 



(B8) 



Consider now the calculation of R{t) for 17q ^ ^ i2 <c 
D~2/3^ Inserting Eq. ^) into ( pb| ) one obtains 



R{t) = 1/2 + N{t)/2 + TTDt{Ii - 2I2) 
Here N{t) is given by Eqs. (||) and (pi). 



(B9) 



, sm X — X cos X . fx\ , 
^1 = / 7^5 / [-j: ) dx 







X^ 



/ smx — xcosx , „, , s 

~ /(O) / 2 dx = /(O), (BIO) 

Jo X 

where the approximate equality holds since t 3> 17^^ and 
the latter equality results from Eq. 3.784.4 in Ref. [ p4| , 
and 

h = r ^^^/ (I) dx^j,+ .h + J3. (Bii) 







X'' 



In Eq. ( p31l| ) I2 is splitted into three integrals, which are 
calculated for t '3> 17g"^ as follows, 



Ji = 



X — sm X „ / X 



X — sm a; 



dx 



= /(0)[7 + sinl-l-Ci(l)] 



(B12) 



J2 - 
and 



sma; „ /x 



/(|)da;«/(G)[Ci(l)-sinl], (B13) 



J3 = 



x\ /(O)^ / 1 



/(0)(lnl7o<-7/2). 



2 ^Hl72i2 



(B14) 



Here i?i() is the integral exponential function a nd C i() is 
t he in tegral cosine The integrals in Eqs. (B12) and 
(B13) were calculated with the help of the Matheniatica 
software p8| , whereas the approximate equality in Eq. 
(B14) results from Ref. p7| , Eq. 5.1.11]. Combining Eqs. 
(p9^-(pl^) yields Eq. 



APPENDIX C: SOLUTION OF EQ. (BIH) 



Inserting Eq. ( |6q ) into Eq. ( |61bD and taking into ac- 
count Eq. (25b), one obtains the equation 



hX 



2„" 



Dh^g 



Dh^ 



vX g' 



2Dh'^ 
^9, 



(CI) 



where X = h{t)Q and the prime denotes the derivative 
of g{X) with respect to X. By the assumption, g{X) 
depends on the time only through X, which means that 
the coefficients in Eq. (CI) should not depend exphcitly 
on the time. This is indeed so, if 



h + uh = -Dh^ 



(C2) 



(C3) 



Equations ( |61b| ) and ( |60| ) and the fact that Lof{^) = 
result in Ko{n 00) /(17). The initial and boundary 



Then Eq. (|Cl| ) becomes 

g" + (l/X + X)g' - {2lX^)g = 0. 
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condition s (pO| ) and ( |65| ) and the latter relation for the [11 
function (|66|)yield the following conditions, [12 

(a)5(0)-0, (h) g{h{0)n) ^ 1, (c) 5(00) = 1. (C4) [13; 



The comparison of conditions (|C4|b) and (|C4p ) shows 
that h{0) = 00. 



Solving Eq. (C2) with the latter condition yields 



(C5) 



The solution of Eq. ( |C3| ) is obtained by the method in 
Ref. l39||. The substitution g{X) = {-Y)''gi{Y), where 



k = I/V2 and Y = -X'^/2, reduces Eq. to the 

degenerate hypergeometric equation |39, Sec. 2.113] 



Yg'l + {2k + l- Y)g[ - kgi 



0, 



(C6) 



where the prime denote s th e derivative with respect to 
Y. The solution of Eq. @ yields 



5(X)=C2C'M(fc,l + 2fc,-C) 



-0, 



(C7) 



where C = -^^72 and C2, C2 are constants. The boundary 
condition (|C^ ) yields C2 — 0, whereas the boundary 
condition (|C4|c) and the asymptotic formula for M() ||2^, 
Eq. 13.5.1] yield Eq. (|6|). As a result, Eqs. @ and 
@ yield Eq. (jeTt). 
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